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ABSTRACT 


Useful  formulal i ons  ol  th('  macroscopic  conservation 
relations  for  a  Ihrcu'-component  mixture  composed  of 
ions,  electrons,  and  neutral  particles  are  developed 
from  the  basic,  forms  of  tlu'  equations  lor  the  individ¬ 
ual  species.  Dk'  di'velopnumt  represents  primarily 
an  orderly  coJT^pil ation  of  liic  work  of  previous  inves¬ 
tigators,  although  th<’  equations  developed  are  some¬ 
what  more  general  than  those  previously  derived.  In 
particular,  all  inertial  and  viscous  terms  are  retained, 
and  the  temperatures  of  the  individual  species  are  not 
required  to  be  equal.  In  addition,  treatment  of  the 
energy  equations  is  included. 
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SECTION  I 


\  ONE- FLUID  FORMULATION 

^^^^The  purpose  of  this  study  is  to  develop,  in  a  systematic  and  general 

way,  useful  formulations  of  the  macroscopic  conservation  relations  (i.  e.  , 

A 

mass,  momentum,  and  energy)  applicable,  in  the  absence  of  radiative  effects, 
to  a  gas  mixture  composed  of  ions,  electrons,  and  neutral  particles,  ^fencer* 
^he  material ;^ofttaincd  in  thio  constitutes  primarily  an  orderly  com¬ 
pilation  of  the  work  of  previous  for  example-  5^hlut^r^^ 

2.3  4.  ^ 

Cowling, - Delcroix,  and  Tiphnerty.^  The  strives 

to  maintain  some  degree  of  generality  so  that  the  various  terms  neglected  in 

other  treatments  are  readily  apparent.  In  particular,  all  inertial  and  viscous 

terms  are  retained,  and  the  temperatures  of  the  species  are  not  required  to 

be  equal.  In  addition,  the  energy  equations,  not  treated  in  previous  works, 

are  included. 

In  Section  II,  the  basic  equations  for  each  species  are  written,  and 
arguments  are  advanced  for  treating  the  momentum  properties  of  the  mixture 
either  as  a  single  fluid  or  a  two-fluid  mixture  and  the  energy  properties  as 
a  one-,  two-,  or  three-temperature  mixture.  Forms  of  the  individual  species 
equations  convenient  for  developing  the  subsequent  formulations  are  presented. 
In  Section  III,  the  one-fluid  formulation  is  developed,  and  in  Section  IV,  the 
two-fluid  formulation  is  derived. 


\ 


1 


SECTION  II 


BASIC  MULTICOMPONENT  EQUATIONS 


A.  Equations  for  Each  Species 

The  governing  equations  for  the  flow  of  a  multicomponent  mixture  com- 
posed  of  species  a  (a  =  1 ,  .  .  .  ,  N)  are  derived  by  taking  moments  of  the 
Boltzmann  equation  appropriate  to  each  species.  In  the  absence  of  any  ex¬ 
ternal  forces  except  those  of  electromagnetic  origin,  the  Boltzmann  equation 
for  species  a  is  written  as 


af  e 

-^+V.Vf+-^(E+V  X  B).V  f 
at  — a  —  m  —  — a  ~  va 


(1) 


where  represents  the  effects  of  interactions  of  members  of  species 

a  with  all  species  presenti  including  other  members  of  species  a. 


Following  the  usual  procedure,  we  obtain  the  transport  equations  (i.  e.  , 

continuity,  momentum,  energy)  for  each  species  by  multiplying  Eq.  (1)  suc- 

2 

cessively  by  m^,  and  1/2  a  integrating  over  velocity 

space,  making  liberal  use  of  the  definition  of  the  mean  value  of  any  function 
^  +  oo 

of  velocity,  =  //;  f^Qd  V^d  V^d  V^.  This  procedure  yields,  after  suit- 

-  oo 

able  rearrangement; 


+  V-(p  V  )  =  io 
at  a— a  a 


(2) 


D“v 

,  _ =5 

a  Dt 


N 


-  Vp  +  V*i  +  Q  (E  +  V  X  B)  +  5]  P  .  -  (A)  V  (3) 

a  -a  a  —  —a  —  j_2  a— a 
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*  f 


Other  useful  forms  of  the  energy  equation  are  obtained  in  terms  of  the  en¬ 
thalpy  and  the  total  enthalpy 


and 


N 


+  z 

j=i 


E  . 
aj 


(4b) 


0°^ 
Pa  T5r 


+  E-  Q  V 
—  a— a 


c 

a 


N 

I 


j=i 


E  . 
QJ 


(4c) 


It  is  to  be  noted  that  can  also  be  expressed  as  3kT^/2m^.  The  physical 
significance  of  the  various  terms  in  the  above  equations  is  well  known  and 
will  not  be  considered  here. 


To  complete  the  formulation,  Maxwell's  equations  are  required: 


8B 

Vx  E  =  - 


(5) 


dE 

^xB  =  pJ  +  p,€ 

—  '^o—  '^o  o  9t 


(6) 
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V-E 


(7) 


=  -Q 

c 

o 

V-B  =  0  (8) 

Also  requirt'd  are  the  equations  of  state  for  each  of  the  species  a: 

It  is  worth  remarking  at  this  point  that  the  basic  unknowns  of  the  prob¬ 
lem  are  p  ,  v  ,  p  ,  E,  B,  u  ,  c  ,  x  »  q  ,  P  K  and  a>  the  entire  set 

a  —a  *ci  —  —  a  q-q  — aj  aj  a 

of  which  comprises  considerably  more  unknowns  than  equations.  We  shall 

therefore  implicitly  assume  that  solutions  of  the  Boltzmann  equation  will 

permit  t  ,  q  ,  P  E  and  (*>  to  be  expressed  in  terms  of  the  other  un- 
^  -a  ^  —aj  aj  a 

knowns  (as  is  the  case  for  dilute  neutral  gases).  The  remaining  set  of  un¬ 
knowns  number  7N  +  6,  two  less  than  the  number  of  equations  available. 

This  is  due  to  the  well-known  fact  that  Eq.  (8)  is  merely  an  initial  condition 
when  used  with  Eq.  (5),  and  the  fact  that  Eqs.  (2),  (6),  and  (7)  are  not  inde¬ 
pendent*  since  (6)  and  (7),  or  (2),  can  be  used  to  obtain  the  conservation  of 
charge  relation.  Hence,  either  Eq.  (5)  or  (8)  and  one  of  Eqs.  (2),  (6),  or  (7) 
may  be  taken  as  superfluous  equations. 

Henceforth,  we  shall  restrict  out  attention  to  three-component  mixtures 
composed  of  ions,  electrons,  and  neutral  particles. 

B.  Useful  Formulations  of  the  Three- Fluid  Equations 

The  equations  for  the  individual  species  listed  in  the  previous  section 
(which  we  shall  refer  to  as  the  three-fluid,  three-temperature  formulation) 
are  generally  not  the  most  convenient  forms  for  dealing  with  specific  three- 
component  flow  problems,  due  primarily  to  certain  simplifications  which 

^It  is  noted  that  the  variables  Q  ,  Q>  and  n  are  defined  in  terms  of 
the  unknowns  listed  here. 
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arise  from  consideration  of  the  coupling  between  the  three  species.  This 
coupling,  as  can  be  seen  from  the  equations  listed,  takes  three  distinct  forms: 
first,  momentum  coupling  due  to  collisions  between  particles  of  different 
species  (the  terms  second,  energy  coupling  due  to  collisions  between 

particles  of  different  species  (the  terms  E  .);  and  third,  electromagnetic 
coupling  due  to  collective  interactions  of  the  charged  species  [the  terms 
Q  (E  +  V  X  B)]  . 

a  —  — Q  — 

If  the  momentum  and  electromagnetic  coupling  between  ions  and  electrons 
is  very  strong,  then  the  number  density  of  each  species  tends  to  be  the  same 
value,  and  the  difference  in  mean  velocities  v^  and  produced  by  mechani¬ 
cal  effects  tends  to  be  small  compared  to  either  or  C^.  This  suggests 
that  it  may  be  useful  to  describe  the  momentum  properties  of  the  ions  and 
electrons  as  a  single  fluid  (the  plasma)  of  mass  density  Pp  ~  P^ 

mean  velocity  p  v  =  p.v.  +  p  v  .  The  mean  velocities  of  the  ions  and  elec- 
^  p  ^1—1  '^e— e 


trons  are  then  described  by  the  diffusional  velocities  w.  and  w  where 

y  ^  -ip  --ep 

.2^  p  w  =0. 

1,  e  ap 

Similarly,  if  the  momentum  coupling  between  ions  and  neutrals  and/or 
electrons  and  neutrals  is  very  strong,  the  differences  in  mean  velocities 


V.  -  v  and  v 
—1  —n  — € 


-  v^  is  small  compared  to  C^,  C^,  or  This  suggests 


that  it  may  be  useful  to  describe  the  momentum  properties  of  ions,  electrons, 

and  neutrals  as  a  single  fluid  of  mass  density  P  “  P^  ■*"  P^  Pj^  mean 

velocity  pv  =  p.v.  +  p  v  +  p  v  .  The  mean  velocities  of  the  individual 
^  “  '^1—1  e— e  ^n— n 

species  are  then  described  by  the  diffusional  velocities  w.,  w  ,  and  w 

i  e  n 

where  .  p  w  =0.  For  this  description  it  is  desired  to  recast  the  in- 

1,  e,  n  a  ^ 

dividual  momentum  equations  into  the  form  of  one  momentum  equation  for  the 

total  mixture,  and  two  additional  diffusion  equations  for  determining  the 

diffusional  velocities.  This  is  the  standard  fluid  mechanical  treatment  of 

5 

multicomponent  mixtures  (see  Hirschfelder,  Curtis,  and  Bird  ). 


If,  however,  the  momentum  coupling  between  the  charged  particles  and 
neutrals  is  not  strong,  it  may  be  desirable  to  describe  the  momentum  prop¬ 
erties  of  the  mixture,  composed  of  (a)  a  plasma  of  mass  density  and  mean 
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velocity  and  (b)  a  neutral  gas  of  mass  density  and  mean  velocity  v^. 
The  mean  velocities  of  ions  and  electrons  are  described  by  the  diffusional 
velocities  w.  and  w  For  this  description  it  is  desired  to  recast  the 
individual  momentum  equations  for  ions  and  electrons  into  one  momentum 
equation  for  the  plasma,  and  one  additional  diffusion  relation  (a  generalized 
Ohm*s  law)  for  determining  the  diffusional  velocities. 


Similar  considerations  may  be  applied  to  the  energy  coupling,  although 

no  exact  reformulations  appear  useful.  However,  if  the  coupling  between 

species  is  sufficiently  strong,  the  assumption  T.  =  T  =  T  =  T  may  be 

1  e  n 

permissible,  and  the  energy  properties  of  the  mixture  may  be  described  by 
a  single  temperature.  It  is  then  desirable  to  combine  the  individual  energy 
equations  into  one  equation  by  eliminating  the  coupling  terms.  Likewise,  if 
the  coupling  between  ions  and  neutrals  is  sufficiently  strong  to  permit  the 
assumption  while  the  electrons  are  weakly  coupled,  then  the 

energy  properties  of  the  mixture  may  be  described  by  two  temperatures, 
and  T^.  In  this  case  it  is  desirable  to  combine  the  energy  equations  for  ions 
and  neutrals  by  eliminating  the  mutual  coupling  term. 


There  exist,  of  course,  other  extremes  in  both  momentum  and  energy 
coupling  between  species  which,  in  theory,  would  suggest  different  formula¬ 
tions.  These  are  not  considered  important  here,  as  they  are  not  physically 
probable. 


To  summarize,  the  momentum  and  energy  properties  of  a  three- 
component  mixture  may  usefially  be  described  as  indicated  below. 


Momentum  Energy 


Three-fluid: 

Pi'  Pe’ 

Pn’  ^i'  :5:n’ 

Three -temperature: 

Two -fluid: 

P  1  V  , 

P  -P 

Pn’  ^ 

Two-  tempe  rature : 

One -fluid: 

p,  V,  w 

e’  ^n’  ^ 

One -temperature: 

T 

C.  Useful  Forms  for  the  Individual  Component  Equations 

As  an  aid  in  recasting  the  equations  in  the  forms  suitable  for  the  descrip¬ 
tions  of  the  previous  section,  the  individual  component  continuity,  momentum. 
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and  energy  equations  will  be  revised  in  the  following  ways:  first,  equations 

for  the  ions  and  electrons  will  be  written  in  terms  of  v  and  the  diffusional 

P 

velocities  w^^  and  and,  second,  equations  for  all  species  will  be  written 

in  terms  of  v  and  the  diffusional  velocities  w.,  w  ,  and  w  . 

—  —1  — e  — n 

In  terms  of  the  first  variables,  the  equations  become  (where  a  =  i  or  e) 


ap 


+V'fpw  ^+V’fpv  \  =  a) 
at  ^  a— p;  a 


(11) 


-  VP 

a  '  —  ap  '  p  '^ap 


a-^p'  *  a—ap  '  — p 

+  V‘i  +Q  (E  +  v  xB)  +  Qw  xB+yp.-ojv 
-an  ^a  -  — p  a— ap  ^  —  aj  a— p 


dP 


Pp  W  [^  (  V  +  *a)] 


-P 

D^v 


(12) 


Pa^ap'^Dt^  "  ■  ’^■^ap  '  Pap 
+  <Iap=^p)  +  +  Xp  ^  B)  -  v  ’  Z  P 


P  j 


2 

+  Z  E  .  +  (a) 

.  aj  a  2 


— ^  (o).  +  Oj)(u  +  €) 
p  '  i  e  a  a 

p 


(13) 


Additional  forms  of  the  energy  equation  are 


Pp  T5r 


(v '  ’•)] 


+€)|+pw  •  =  -  V  •  £  +  V  •  Vp 

ap  a/  I  ^a— ap  Dt  ^ap  — p  *^ap 


+  (l  )  +  Qw  -(E  +  v  xB)-v*5]P. 

-ap  ^p'  a-ap  -p  -  -p  r -aj 


+  Ze.+-^  +  w  -§■ 
j  aj  at  a  2 


— ^  (o).  +  U)  ) 
p  i  e 

P 


^  ap  a 

(14) 
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and 


dP 
Pp  w 


+  c  +  -# 

a  Z 


4-  ^ 

pp  :Dr 


—  (w  •  V  ) 

P„  -ap  -P 


+  V-\  \  p  v^w  )=-V"q  +V*(v-T  ) 

\Z  *^0  l>~ci/  — ip  ~P  ®P 


+  E-Q  V  +  E.Q  W_  +  I  +  -^ 


—  a— p  —  a— ap 


“  ~  (ci).  +  a>  )  I  (w  •  V  )  +  ( +  h  «  1 1  /I  c\ 

p  '  1  e'  -ap  -p'  12  ap  all  (15) 

A  ^ 

It  is  to  be  noted  that  in  the  above  equations  the  quantities  P^p»  '^ap’^ap'  -ap’ 
and  are  defined  in  terms  of  tl;ie  appropriate  mean  values  of  the  difference 

in  particle  velocity  and  mean  plasma  velocity  (C^p),  instead  of  the  difference 
in  particle  velocity  and  mean  species  velocity  (C^). 

In  terms  of  the  second  variables,  the  equations  become  (where 
a  =  i,  e,  or  n) 


+  V-  (p^w^)  +  V.  (p^v)  =  (16) 

D 

Pa  T5t  Dt  <PQ^iS^a>  Pa^a  Pa  ^  = 

-  Vp  ^+V'i  +Q  (E  +  vxB)  +  Qw  xB  +  5]P.-a)V 

•^at  -at  a  —  —  —  a— a  —  ^ 

(17) 
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D  r^a  /A  vl  .  Dv 

I  —  (u.+€  )|  +  pw  --pir 
Dt  I  p  at  a  I  a— a  Dt 


-  -  Pat 


+  Qw"(E  +  vxB)  +  Ze.  -  vZ  P  .  +  CO  ^  (18) 

^Q-a  -  -  j  ~  j  ^  ^ 


Alternative  forms  of  the  energy  equation  are 

Dv 


D  /L  1  \  I  Dv  _  J.  r-r 

P  ^  [T  ^^at  +  'a]  +  Pa^a- dT  =  ' 


and 


+  (l„.:Sy.)  +  Q„w  •  (E  +  V  X  B)  +  z  E  . 
-at  —  a — a  —  ^  aj 


V  Z  P  . 
-  j 


2  ap  . 

L  V  '^at 


(19) 


Dt  [  p 


a)]  ^  P  m  [t  -a  -]  +  ^a] 

-  V-  q  *  +  V*  (v'  T  .)  +  E-  Q  V  +  E-  Q  w 
^at  '—  =at  —  a—  —  a— a 


>  .  +  V  +  * 

L  at  2 


8P  f 

+  I  E  :  +  ^  (^“1 

^  «J  3t 

A  A 

Again  it  is  to  be  noted  that  the  quantities  u^^,  h^^,  and  are 

defined  in  terms  of  the  appropriate  mean  values  of  the  difference  in  particle 
velocity  and  mean  mixture  velocity 
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SECTION  in 


ONE -FLUID  P^ORMULATION 

A.  Basic  Fluid  Mechanical  Formulation 

Considering  the  mixture  as  a  single  fluid  with  momentum  p^roper ties 
described  by  p  and  v,  we  obtain  the  equation  of  continuity  by  ^ldding  the  three 
Eqs.  (16)  and  making  use  of  the  facts  that  ^  ^ 

+  V-  (pv)  =  0  (21) 

The  conservation  of  the  individual  species  may  be  expressed  directly  by  (16) 
where,  it  is  to  be  remembered,  only  one  such  equation  is  required  in  com¬ 
bination  with  Eq.  (21)  and  Maxwell’s  equations. 

The  momentum  equation  for  the  mixture  is  obtained  by  adding  the 
Eqs.(17): 


Dv 

p  ^  -  Vp^  +  +  Q(E  +  V  X  B)  4  X  B  (22) 


A  form  for  the  diffusional  relations  useful  for  comparison  with  standard 
fluid  mechanical  developments  is  obtained  by  eliminating  the  Dv/Dt  term 
from  Eqs.  (17)  and  (22).  (Only  two  of  these  relations  are  independent.  ) 


z 

j 


p  . 


+  ’Pat  + 


Q^j  (E  +  V  X  B) 


+  f—  j  -  Qw)xB+a)v  +  —  /  y.  ^  ) 
\p  ^  a-a/  -  a-  p  '  =V 


(V- 


lat) 


(23) 
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These  equiitions  merely  represent  the  conservation  of  the  difference  in 
momentum  of  a  single  species  and  the  mixture;  the  left-hand  side  is  the 
friction  force  acting  on  a  single  species  due  to  interaction  with  the  other 
species,  and  the  right-hand  side  represents  the  summation  of  effects  which 
require  this  interaction  for  dynamical  equilibrium.  In  addition,  these  equa¬ 
tions  are  of  the  form  obtained  by  Hirschfelder ,  Curtiss,  and  Bird^  by  the 
standard  perturbation  solution  of  the  Boltzmann  equations  for  the  species, 
in  which  it  is  assumed  that: 

1.  All  species  are  of  equal  temperature; 

2.  The  velocity  distribution  function  of  each  species  is  perturbed 
slightly  from  the  Maxwellian  one  corresponding  to  the  mean 
properties  of  the  mixture; 

3.  The  forces  on  the  particles  arising  due  to  E  and  B  do  not 
directly  affect  the  distribution  function; 

4.  There  is  no  chemical  reaction  (cj  =  0). 

a 

5 

The  precise  equations  obtained  by  Hirschfelder,  Curtiss,  and  Bird, 
when  written  in  the  form  of  Eqs.(23),  are: 


p  .  =  -  ^  Vp  +  Vp  +  (^  Q  -  Q  )  (E  +  V  X  B) 
.  —aj  p  ^t  ^at  \  p  a/  —  —  — 


+  (^  i  -  Q  A  5) 


(24) 


where 


P  .  = 

-aj 


n  n. 

J 


„  D^.(l) 


w  )  + 
—  a 


n^DjU) 


(25) 


U 


The  usefulness  of  comp^lrin^  Eqs.  (Z4)  and  (Z5)  with  Eq.  (23)  is  that  it  makes 
evident  the  effect  of  the  approximations  made  in  obtaining  (24)  and  (25)  [the 
last  five  terms  on  the  right-hand  side  of  (23)  are  neglected]  ,  and  it  further 
yields  the  form  of  the  friction  terms  (subject,  of  course,  to  the  previous 

assumptions).  It  is  pointed  out  that  the  terms  neglected  in  Eq.  (23)  are  a 
direct  consequence  of  the  .assumption  that  the  velocity  distribution  function  of 
each  species  suffers  only  small  departures  from  the  Maxwellian  one  corre¬ 
sponding  to  the  mean  properties  of  the  mixture.  If  the  coupling  between 
species  is  strong,  we  may  expect  the  mean  properties  of  each  species  to 
differ  little  from  the  mean  mixture  properties;  this,  of  course,  represents 
the  situation  in  which  the  one -fluid  model  is  convenient. 

For  the  one -temperature  model,  the  energy  equations  are  obtained  by 
adding  the  Eqs.  (18),  (19)»  or  (20),  respectively: 


D  fv  P 

m  [I  -  V 


.  +  € 
at  a  i 


V  qt  -  Pj(V  v)  + 


+  (E  +  V  X  B) 


or 


p  m  (\t  +  *a)j  =  -  +  rzPt  + 

+  (E  +  V  X  B)  + 


or 

P  m  [2^^  (^at  +  *a  +  v)]  =  - 

9Pt 

+  E.Qv  + 


(27) 
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Similarly,  the  appropriate  energy  equations  for  the  two -temperature 
model  are  obtained  by  adding  Eqs.  (18),  (19)i  or  (20)  for  ions  and  neutrals, 
respectively: 


P  Dt 


i,n 


Dv 

Pe^  DF  =  - 


it  ' 


-  ^Pit  +  Pnt^  ^  [<Iit  + 

+  Q.w.  (E  +  V  X  B)  +  E.  +  E  -  v(P.  +  P  ) 

1—1  '  le  ne  —  — le  — ne' 

2 

+  (“i  +  ^  (29) 


or 


D 

P  Dt 


^  T  (^at  +  ‘  a) 


i,n 


Pe^  m  =  ■  ^  V  ^  rv(P4,  +  P„,) 


It  ■  ^nt' 


+  Rt  *  +  T  J:Vvl  +  Q.w.  (E  +  V  X  B)  +  E.  +  E 
-it  =nt'  =J  1—1  le  ne 


a 

v(P.  +  P  )  +  (cO.  +  O)  )  +  *5^  (p,.  +  P  x)  (30) 

—  le  — ne'  '  i  n'  2  3t  '^it  ^nt'  '  ' 
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( 


or 


D 

P  Dt 


i,n 


h  +  c  '\-  -rr 
at  a  Z 


-  +  Ent)  +  +  Int)] 


+  E-  Q  V  +  E-  Q.w.  +  E.  +  E  + 
—  o:-  —  1—1  le  ne 


^<Pit  +  Pnt^ 
dt 


(31) 


The  energy  equation  for  the  electroUvS  is  obtained  directly  from  Eqs.  (18), 
(19),  or  (20). 

For  the  three -temperature  modtd,  the  appropriate  energy  equations 
are  obtained  directly  from  Eqs.  (18),  (19),  or  (20). 

B.  Form  of  the  Friction  Forces,  P  . 

_ 1152 


It  is  not  profitable  to  proceed  with  the  one -fluid  model  without  intro¬ 
ducing  the  form  of  the  momentum  coupling  terms  P  ..  The  form  adopted 
here,  for  reasons  to  be  discussed,  is  as  follows: 


p. 

=  -  P  . 

Tl .  •  (w 
die  e 

-  w.)  +  P. 

(32) 

—  le 

—  ei 

— i'  — le 

P. 

=  -  P  .  = 

Pin 

-  w.)  +  P. 

(33) 

—  in 

—  ni 

—  1  —  in 

P 

=  -  P 

Pen 

-  w  )  +  P 
— n  — ne 

(34) 

—  ne 

—  en 

^ie 

and  T)  are  of  the  form 

Jen 

'H 

.  =  T1  .  I 

"  i 

(35) 

aj  'laj  = 

IPJ 

'laj  — 
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where  X  i®  unit  tensor  and  n  is  the  unit  vector  in  the  direction  of  the 
magnetic  field. 


The  form  of  these  relations  suggested  by  comparison  with  the  standard 

T 

fluid  mechanical  treatment  is,  from  Eq.  (Z5),  P  .  =  n  .(w.  -  w  )  +  n  .VT. 

Qj  -a'  •aj"' 

It  is  recalled,  however,  that  this  form  includes  the  assumption  that  electric 

and  magnetic  fields  do  not  directly  affect  the  velocity  distribution  function; 

this  has  been  shown  by  other  investigators  (e.  g.  ,  Chapman  and  Cowling, 

7  8 

Marshall,  and  Braginski  )  to  be  imprecise.  In  particular,  Braginski 
obtains,  for  a  mixture  of  ions  and  electrons,  the  following  form: 


P  .  =  n  .*(w.  -  w  ) 
-aj  Jaj  “J  -a' 


VT  +  rf 


1*  .  Mw.  -  w  )  X  —2-1 

Vj  ,3,J 

1^.  /vT  X 

\—  iBlj 


(36) 


where  it  has  been  assumed  that  the  velocity  distribution  function  is  unaffected 
by  electric  fields.  The  form  of  the  tensors  r|  is 

^  I  -  (t1||  -  nn  (37) 

Clearly,  if  the  electric  fields  were  properly  accounted  for,  a  term  propor¬ 
tional  to  E  is  expected  in  Eq.  (36). 

Examination  of  Braginskies  results  reveals  that  over  all  ranges  of 
magnetic  field  strength,  in  (37)  varies  from  qj|to  2q||;  is  always 
less  than  20  percent  of  the  value  of  q||.  Furthermore,  for  moderate  magnetic 
field  strengths  1 1»  velocity  distribution  function  is 

only  slightly  affected  by  the  magnetic  field.  These  considerations  suggest 
that  the  forms  of  P^j  in  Eqs.  (32)  through  (34)  give  proper  emphasis  to  the 
most  important  effects. 
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C.  Onc-Kluici  Ko  rniul*!  I  ion  in  'J'('rms  of  p,  f,  Q,  v,  j,  and 

"  1 .  Relationships  Between  Variables 

It  is  found  most  convenient  to  cast  the  one- fluid  formulation  in 
terms  of  the  variables  p,  f,  Q,  v,  j,  and  We  list  here  the  other 

important  variables  in  terms  of  these: 


Pi  =  “2 


n 


a,  fp 


(38a) 

(38b) 


o-  5.-0,  rn  e. 

a/  -  (1  -  f)p 

Z  6  m.e  ' 

1  e 


(38c) 


Q. 


❖ 


a 


2 


e. 

1 


m, 

1 


(1 


f)p 


(38d) 


❖ 

Q  =  Q  -  Q,  —  (1  -  f)p 
e  2  m. 

1 


w. 
—  1 


i  (1  -  fa^) 

*  -^7  (1  -  f)p  ~y~  (1  -  f)  -A 


Q-  m 
5  e 


e 


w  = 
—  n 


ie  aSe 

04  (1  -  ^^2^ 

a^f 


(38e) 


(38f) 


(38g) 


w 
—  e 


asi 


* 

a^a.  (1 


Sjt  ❖❖ 
^2  % 


e. 

—  (1 

m.  ' 

1 


-  f)p 


(1  -  f)  a 


fa*) 

V  Xa 


(38h) 
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where 


Furthermore 


a 


2 


1  + 


m 
_ e 

m. 

1 


f 


m  e . 
e  1 


m.e 
1  e 


(1 


f) 


(39a) 


(39b) 


(39c) 


(39d) 


(39e) 

(39f) 


,  we  may  write 


w.  -  w  = 
—  1  — e 


A. 

le 


m. 
1 


^  e.  (1  -  f)p  j-  ^  °ie 


“A  / 


(40a) 


w  -  w 
— n  — e 


Qc  m. 

5  1  . 

“5  ^  e.  (1  -  f)p  i- 

Q,  1  ^ 

L  D 


) 

ne 

f  -A 


(40b) 
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w 
—  n 


w, 
—  1 


(40c ) 


where 


m 


L 


m.e 
1  e 


(4la) 


D. 


le 


5  4 


a,, 

Li 


(1  -  fa^) 

(1  -  r) 


D 


) 

ne 


T 


t  ^  +  j  “IT"^  .  f) 


a- 


D. 

in 


ci^a^  ^5^4  ~ 

—  -  f  -1^  +  f  ^ - (T  -  f)  ■ 

^3  ^3  ^2 


(41b) 


(41c) 


(41d) 


For  the  sake  of  brevity,  we  shall  now  make  the  standard  assumption  of 

4 

charge  neutrality  (see  Lehnert  )  and  further  assume  that  m^/m^  <  <  1  and 

|m^e^/m^e^[  <  <  1.  In  the  absence  of  the  assumption  of  charge  neutrality, 

there  appears  to  be  no  advantage  in  the  one -fluid  formulation.  With  these 

assumptions,  the  coefficients  a,  A,  and  D  in  the  above  equations  become 

unity.  ^  We  shall  neglect  the  unspecified  terms  in  the  friction  forces  between 

species  (P  .). 

-cij 


^In  order  to  obtain  the  correct  limits,  it  is  necessary  to  perform  the 
subsequent  algebra  by  retaining  the  coefficients  a,  A,  and  D,  and  then 
invoking  the  previous  assumptions. 
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Z,  Generalized  Ohm *8  Law  and  Ambipolar  Diffusion  Equation 


In  terms  of  the  new  variables,  it  is  convenient  first  to  extract  an 
Ohm's  law  for  j_  and  a  diffusion  equation  for  from  the  three  dependent 
relations  (23),  where  Eqs.  (32)  through  (34)  with  (40a-c)  are  used  to  eliminate 
P^..  To  obtain  a  diffusion  equation  for  v^,  Eq.  (23)  for  neutrals  is  written 
in  terms  of  the  new  variables,  with  the  result  that 


(^en||  I  ^in) 


V.  =  -  f  Vp.  ,  +  Vp  .  +  f(j  X  B) 

_A||  j_  -^t||  1  j_  J-  -  ,|j_ 

e.  \  m. 

^  ^inj  ej(l  -  f)p  j-||  J_ 


+  h 


+  CO 


V  + 

-II.  1 


T.) 

=t'll 


.1 


(V 


X 


f) 


(42) 


To  obtain  a  generalized  Ohm's  law,  Eq.  (23)  for  ions  is  subtracted  from  the 
similar  one  for  electrons,  and  Eq.  (42)  is  employed  to  eliminate  v^.  The 
result  is 
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4,1  . 
'^11,1 


m. 

=  (E  +  V  X  D)„  j_  -  e,  p(i  ■  ly  (*  -  2  fV||,j_  )  (I  ^)|1,1 


e  1 


II 1 


^  *-pVx  "  ^  ['“  ■  'll.x’  '  ''I'.J-  “*]  'i  ■  *'*  -I'-L 

mm.  f  Di  *1 

jro  [pt'  -  ^Xa]  +  PO  ■  '>  [xa<’  x'  "  '-a  ’>x]|||  ^ 

Jp  +  It  +  (1  +  -T  ^  -  T  jm  ■  ^^-a] 


p(l  -  f)  +  {lA-^)v]jj  ^  I 


(43) 


where 


'eiii 


^11,1 


(44) 


and 


(”11,1)''  '[^7^^] 


( 


■  V 


(nn  e.  \ 


(45a) 
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nu 


-p'll.l  ‘‘i  po  -  '■) 


^11,1  ■  ^11,1 


e  .m 
i  e 


^  j-i;r  •  <■  -  (yE,,Y_^ 


e  1 


(45b) 


""^’ll,X  ‘  '  i,  pTi’-“i)  [''ll.x*'  ■  Jt’ii.i-  i'ii,x'’  Jit>i|,i 


e.m 

+  (V  t  )  +  (1  -  Y  )  (V  t 

®e""i  ».-L 


F  =  i-4—  Vp,  -  -1  Vp.,  -  I  Vp 


t^ia] 


(45c) 


(45d) 


-p  f  f  ’nt  f  "^et 

r.  =  -  (^)  ('^•1,1  +  T  +  T  <''■5,.' 


Schluter^  has  derived  equations  similar  to  Eq.  (43),  in  which  and 

^ie  are  scalars  and  the  viscous  stresses  are  neglected.  Cowling  has 
employed  Schluter's  assumptions  in  deriving  similar  equations,  assuming 
in  addition  that  the  temperatures  of  all  components  are  equal  and  that  the 
inertial  terms  involving  j_  and  v^  are  neglected.  Briefly,  the  terms  in  this 
equation  represent  current-producing  mechanisms  due  to:  (a)  electric  field, 
(b)  Hall  effect,  (c)  different  motions  of  ions  and  electrons  caused  by  pressure 
gradients,  (d)  different  motions  of  ions  and  electrons  caused  by  viscous 
stresses,  (e)  generation  of  ions  and  electrons  of  unequal  velocities,  (f)  dif¬ 
ferent  inertial  forces  acting  on  ions  and  electrons,  (g)  different  inertial 
forces  acting  in  the  ion-electron  mixture  and  the  neutral  particles  and 
(h)  relative  motion  between  ions  and  neutral  particles.  The  form  of  this 
equation  for  a  fully  ionized  gas  is  readily  obtained  by  setting 

f  =  v.  =a).  =  0)  =0. 

—A  1  e 
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In  Eq.  (42),  the  various  terms  represent  motion  of  the  plasma  rela¬ 
tive  to  the  mixture  due  to:  (a)  action  of  pressure  gradients  (this  includes 
both  mass  diffusion  and  pressure  diffusion),  (b)  the  Lorentz  force  acting  on 
the  plasma,  (c)  currents  flowing  in  the  plasma,  (d)  generation  of  neutral 
species  (i.e.  ,  recombination  of  ions  and  electrons),  (e)  action  of  viscous 
stresses,  and  (f)  inertial  effects.  Schluter  has  also  derived  a  similar 
equation,  and  has  pointed  out  that,  if  the  last  five  effects  are  absent,  the 
temperatures  of  all  components  are  equal,  and  the  gas  is  slightly  ionized 
(f  <  <  1),  then  Eq.  (42)  reduces  to  the  usual  ambipolar  diffusion  equation: 


^A. 


,1 


2kT  \  „ 

- 1 - 

’en||  ^  ’^in/ 


3.  Continuity  and  Momentum  Equations 

Incorporation  of  the  previous  assumptions  leaves  the  continuity 
equation  unchanged;  hence,  from  Eq,  (21) 


1^  +  V  (pv)  =  0 


(46) 


Recalling  that  only  one  species-conservation  relation  is  necessary,  we 
obtain  from  Eq.  (16),  for  the  neutral  species: 

^  (fp)  +  V-  j^p(l  -  "  “n  ■ 

The  momentum  equation  is  obtained  from  Eq.  (22)  with  Q  =  0 

Dv 

P  ^  ^  +  (i  X  B)  (48) 
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4.  lOnc^rcyy  Equations 

1'hc  energy  equations  for  the  three -temperature  model  are 
obtained  directly  from  Eqs.  (18),  (19)»  or  (ZO)  using  the  relations  (38) 
through  (41).  The  form  corresponding  to  Eq.  (ZO)  becomes,  for  ions: 


Op. 


+  E.  +  E.  + 
—  in  “  le 


it 


at 


(49a) 


for  electrons: 
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and»  for  neutrals: 


Relations  corresponding  to  Eqs.  (18)  and  (19)  can  similarly  be  obtained. 

The  energy  equations  for  the  two-temperature  model  arc  obtained 
from  Eqs.  (29)i  (30),  or  (31).  For  the  mixture  of  ions  and  neutrals,  the 
form  corresponding  to  Eq.  (31)  becomes 


+  E,  +  E  + 
—  ie  — ne 


»'Pit  ^  Pnt' 


The  energy  equation  for  electrons  for  this  case  is  again  (49b). 
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The  one -temperature  formulation  for  the  energy  equation  is  obtained 
from  Eqs.  (26),  (27),  or  (28).  The  form  corresponding  to  Eq.  (28)  becomes 

f  m[<>  -  «  'Nt  +  s'  -  (^t  + ',)  +  ‘(Kt  +  <n)  +  x]  ^ 


-  V-q^  +  V-(v.  X^)  +  + 


^Pt 

at 


(51) 


It  is  pointed  out  that  the  relationships  of  the  energy  flux  vectors  q^^, 

and  to  the  corresponding  vectors  relative  to  the  individual  species  q^ 
and  3^^  are  given  by: 


Sit  =ai 


m 

"1  r  m  ‘ 

/A  \ 

e  ; 

lit  +  P(1  -  f)vA  -  I 

J  L  ®  . 

{\t  +  *i) 

-A  e^(l  -  f)p 

^t 


■^e  1^— A  eTO  ^  -et 

Sn  +  ^A  ?nt  -  -  ^)^A  (^nt  +  'n) 


h  4  €  ) 
et  e/ 


where  the  terms  represent  (a)  energy  flux  due  to  transport  within  species, 

(b)  work  performed  by  shear  forces  due  to  relative  motions  of  species,  and 

(c)  transport  of  energy  due  to  relative  motion  of  species.  Of  these,  the 
term 


e 

is  the  most  widely  recognized. 
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5.  Concluding  Remarks  on  the  One -Fluid  Formulation 

In  addition  to  the  continuity,  momentum,  energy,  and  diffusion 
equations  presented  previously,  Maxwell's  Eqs.  (5)  through  (7) 

[Eq.  (7)  is  unnecessary  if  Q  =  O]  and  the  equations  of  state  (9)  through  (10) 
are  required.  In  the  latter  connection,  it  should  again  be  noted  that  the 
quantities  p^,  p^^,  p^^,  p^^,  et  cetera,  differ  from  the  corresponding  ones 
^i*  ^e*  ^n*  cetera,  by  factors  containing  the  diffusional  velocities 
squared  hence  these  factors  must  be  negligible  compared  to  if 

any  confidence  is  to  be  placed  in  the  equations  of  state. 
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SECTION  IV 


TWO-  FLUID  FORMULATION 

A.  Basic'  Formulation 

The  two-fluid  mc^dcl  (i.e.  ,  plasma  and  neutral  gas)  is  best  formulatc^d 
in  terms  of  the?  variables  p  ,  p  ,  Q,  v  ,  v  ,  and  i  •  For  future  referemee, 
the  relations  between  the*  diffusional  veloc  ities  and  plasma  current  arc* 


w 

-ep 


P. 
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and  net  charge  density  are 
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(52a) 


(52b) 


(52c) 


(52d) 
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The  continuity  equation  for  the  plasma  is  obtained  by  adding  Eqs.  (11) 

for  ions  and  electrons  (noting  that  Z  P  “0) 

i,  c  ^ 


dp 


=  V.(p  V  )  =  O). 
at  p  — p  1 


+  0) 


The  continuity  equation  for  the  neutral  gas  is  obtained  directly  from  Eq.  (2) 

9P. 


-5^  +  V*  (p  V  )  =  O) 
0t  n  — n'  n 


(54) 


Addition  of  £qs«  (12)  for  ions  and  for  electrons  yields  the  momentum 
equation  for  the  plasma 

D^v 

Pp  -  VPp  +  V-Ip  t  Q(E  +  Vp  X  B)  +  (i^  X  B) 


+  P*  +P  -(a),+u))v 
—  in  — en  i  e  — p 


(55) 


The  momentum  equation  for  the  neutral  gas  is  obtained  directly  from  Eq.  (3) 


n 


d\ 

— n 

T5r 


-  Vp  +  (V-  T  )  -  P. 

'  =n'  —  1 


P  -  CO  V 
in  —  en  n  — n 


(56) 


The  appropriate  form  of  a  generalized  Ohm*s  law  is  obtained  by 
multiplying  Eq.  (12)  for  ions  by  e./m^,  multiplying  Eq.  (12)  for  electrons 
by  e^/m^,  adding  the  two  results »  eliminating  the  term  in  this  result  con¬ 
taining  D^^/Dt  with  Eq.  (55),  and  introducing  the  relations  (52).  This 
results  in 
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-1^  *i.p  f’  V  * 


V‘  T.  ) 
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(  e  \  ^  n^'. 


ep-  \  P.  P^  Pp 


{—  +  —  -  — il  X  B)  +  (—  -  P.  +  (—  -  P 

Pp/  ‘^P  -  Pp/-^"  V”"e  Pp/-^" 

fli  -  P.  +  -9.  (..  +  c  )v 
V^i  Pp  '  ®~P 


As  will  be  seen  subsequently,  this  equation  takes  on  a  more  recognizable 
form  upon  substitution  of  appropriate  forms  of  the  coupling  coefficients 


P.  ,  P  ,  and  P.  . 
—  in  —  en  —  le 


The  energy  equations  for  ions  and  for  electrons  for  the  three- 
temperature  model  are  obtained  directly  from  Eqs.  (13),  (14),  or  (15).  The 
form  corresponding  to  Eq.  (15)  is  as  follows  for  ions: 
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For  electrons: 


+ 


(E.v  ) 
- P 


+  E 

en 


E. 

le 


8p 

+  -  (w.  +  O)^) 


V 

-p 


For  neutrals,  the  energy  equation  comes  from  Eqs.  (4a),  (4b),  or  (4c). 
From  (4c), 


Forms  corresponding  to  Eqs.  (13)  and  (14)  for  ions  and  electrons  and 

Eqs.  (4a)  and  (4b)  for  neutrals  can  similarly  be  obtained.  Relations  (52c-d) 

may  be  used  to  express  the  factors  p./p  >  p  /p  in  terms  of  p  and  Q,  if 

1  p  e  p  p 

desired. 
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For  I  he  two -temperature  modeh  the  energy  equation  for  the  ion -neutral 
mixture  is  obtained  by  merely  adding  Eqs.  (58a)  and  (58c);  the  equation  for 
electrons  remains  (58b). 


For  the  one -temperature  model,  the  energy  equation  for  the  mixturf‘ 
is  obtained  by  adding  Eqs.  (58a).  (58b),  and  (58c) 
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Pn  T5r 


+  ( 


n 


n 


i5r 


Pi  A  Pe  A 

p  '  ip  i'  p  '  ep 
.  P  P 


2 

V 
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(%  +  a„)  +  (''j 


Ip' 


+  V 


(v  •  T  ) 
— n  — n 


“f  E*  Qv 
“  “P 


+  E 


«Pp 

~sr 


dp 

~W 


(w.  + 


w  )  ^ 

e'  p 
_  P 


A 

(h; 


IP 


Pg  A 

—  (h 
Pp  ep 


+  V 


> 


(59) 


Similar  equations  written  in  terms  of  either  enthalpy  or  internal  energy  can 
also  be  obtained. 


B.  Approximate  Two-Fluid  Formulation 

If,  as  in  the  case  of  the  one -fluid  formulation,  we  introduce  the  Eqs. 

(32)  through  (35)  for  the  form  of  the  momentum  coupling  terms,  and  if, 

further,  we  assume  that  m  /m,  «  1,  I  e.  m  /e  m.|  «  1,  and 

^  e  1  I  1  e  e  il 

I  m^  ^pl  equations  in  the  preceding  section  can  be  written 

in  the  following  form. 

Continuity: 


9p„ 

+  V-  (p  V  )  =  0).  +  (a)  (60) 

m  P  “"P  in  '  ' 


^It  is  to  be  noted  that  this  assumption  does  not  require  macroscopic 
charge  neutrality. 
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(61) 


®Pn 

_  “  +  V  (p  V  )  =  w 
-n'  n 


Momentum: 


D*'v 


p 


^  .  Vp  +  V-l  +  Q(E  +  V  X  B)  +  (j_  X  B) 


=P 


-  (CO.  +  co^)v  +  (v^  -  Vp)  +  (v^ 


-  yJ 


+  j  +  6” 

P  e  ^ 
e  ^ 


(*)» 


j  +  p-  +  p  „ 

en  —in  — en 


D  V 
— n 

Pn  “DT- 


=  -  Vp^  +  (V-tJ  -  ‘^n  -  ^in  "  V  '  ?en 


(v  -  V  ) 


Ti.  m 
'in  e 

P  ® 


^  ^  Pp)  '  “ 


.  P 
in  —  en 


(63) 


whe  re 


'••  ■  (■  -  9i)' 


(64) 
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Ohm* 8  law  (after  some  rearrange' mi' nt): 


where 
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Energy  (throe -to  inpc' rntu  re  moch’l): 


+  E. 
in 


+  E. 
le 


(co.  +  O)  ) 
'  1  e' 


(67a) 


34 


n 


D 

T3t 


-  V-q  +  V 


n 


<v  -  T  ) 
— n  =n' 


8p 

~sr 


-w  (fi  +€  +-3]-E  -E.  {67c) 

n  \  n  n  /  en  in  '  ' 

Energy  (one*temperature  model): 


Equations  (60)  through  (66),  in  the  absence  of  viscous  effects  and  chemical 

change,  assuming  charge  neutrality  and  assuming  ^  ,  q.  to  be  scalars, 

have  been  thoroughly  discussed  by  Lehnert.  In  particular,  Lehnert 

assumes  n  =  m  n  v  ,  n.  =  rn.n.v.  ,  and  n.  =  m  n  v.  ,  where 
'en  e  e  en  'm  i  i  in  'le  e  e  le 

the  v's  are  suitably  defined  collision  frequencies. 
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NOMP:NC:j^A'rURE 


a 


D‘^/Dt 


nicignctic  intensity 

peculiar  velocity  of  a  particle  of  species  a;  C  =  V  -  v 
^  ^  ^  ^  --a— a— a 

peculiar  velocity  of  a  particle  of  species  a  relative  to  mass 

velocity  of  plasma:  C  =  V  -  v 

^  ^  —  ap  —  a  —p 

peculiar  velocity  of  a  particle  of  species  a  relative  to  mass 

velocity  of  mixture;  C  .  =  V  -  v 

^  —  at  —  a  — 

binary  diffusion  coefficient  of  species  a  and  j 
thermal  diffusion  coefficient  of  species  a 
substantial  derivative  operator  relative  to  species  a; 

0*^701  =  (8 /at)  +  V  -v 

—a 


D/Dt 


E 


E  . 
aj 


e 


a 


substantial  derivative  operator  relative  to  mixture: 

D/Dt  =  (0/at)  +  v-V 

electric  field 

rate  of  energy  transfer  per  unit  volume  from  species  j  to 
species  a 

charge  of  a  particle  of  species  a 


f 

f 

a 

A 

h 

a 


molar  fraction  of  neutral  particles;  f  =  n^/(n^  +  n^) 

velocity  distribution  function  of  species  a 

specific  perfect  gas  enthalpy  of  species  a  relative  to  mass 

A  A 

velocity  of  species;  h^  =  u^  + 

specific  perfect  gas  enthalpy  of  species  a  relative  to  mass 

A  A  . 

velocity  of  plasma;  h  =  u  +  p  ip 
^  ^  ap  ap  ^ap 
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NOMENCLATURE  (Continued) 


at 


J 

L 


k 

m 

a 

N 

n 

n 

a 

P  . 

“CIJ 


specific  perfef:t  gas  enthalpy  of  species  a  relative  to  mass 

A  A 

velocity  of  mixture;  h  .  •  u  ,  +  p  /p 
^  at  at  ‘  at 

N 

total  current  density;  J  -  T  n  e  v 
^  a  a— a 

a~  1 

conduction  current  density  relative  to  mixture, 

j  “  n.e.w.  +  n  e  w 
11—1  e  e— e 

conduction  current  density  relative  to  plasma, 

j  =  n, e.w .  +  n  e  w 

*^p  1  1— ip  e  e— ep 

Boltzmann's  constant 

mass  of  particle  of  species  a 

number  of  species  present 
total  number  density  of  mixture 
number  density  of  species  a 

rate  of  momentum  transfer  per  unit  volume  from  species  j 
to  species  a 

total  plasma  pressure  p  =  p.  +  p 

N 

total  pressure  of  mixture  p^  ^  Y, 

a=l 

partial  pressure  of  species  a  relative  to  mass  velocity  of 

2 

species;  p  =  1  /3  p  C 
^  '^a  a  a 

partial  pressure  of  species  a  relative  to  mass  velocity  of 

plasma;  p  =  1  /3  p  C^ 

^  '^ap  ^a  ap 
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NOMENCLATURE  (Continued) 


4 


*'at 

Q 

Q 

a 

St 

aa 

^ap 


^at 


T 

T 


partial  pressure  of  species  a  relative  to  mass  velocity  of 

mixture;  p^  =  l/3p  C. 

at  ^  a  at 

total  charge  density  of  mixture;  Q  =  X  Q 

a  ^ 


charge  density  of  species  a;  Q  =  n  e 

a  a  a 

energy  flux  vector  of  mixture; 


energy  flux  vector  of  species  a  relative  to  mass  velocity  of 
2  _ 

species;  q  =  pCC  +  p€C 
^a  a— a— a  a  a—  a 

energy  flux  vector  of  species  a  relative  to  mass  velocity  of  plasma 
2 

pCC  +p€W  =q  -  w  .T 
'^ap  ap  ap  ^a  or-ap  —a  —  ap  =ap 

4*pw  [u  +(p  /p)  +  €] 

a— ap  ^  ap  ^ap  a  a** 

energy  flux  vector  of  species  a  relative  to  mass  velocity  of 
mixture; 

2 

^4.=pC.C.  +  p€W  =  q  -  w.i. 

^at  at  — at  ^a  a— a  ^a  —a  -at 

+  pw  [u.  +  (p  /p)  +  €  ] 

^a— a**  at  ^ap  ^a  a'* 

temperature  of  mixture 

temperature  of  ion-neutral  mixture 


temperature  of  species  a;  =  l/3(m^/k)C^ 
time 

mean  translational  kinetic  energy  per  unit  mass  of  a  particle  of 

species  a  relative  to  mass  velocity  of  species; 

u  =  1/2C^  =  3/2  (k/m  )T 
a  a  a  a 
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NOMENCLATURE  (Continued) 


«  I 


A 

U 


ap 


at 


mean  translational  kinetic  energy  per  unit  mass  of  a  particle  of 

A 

species  a  relative  to  mass  velocity  of  plasma;  u^^  =  ^/^C^p 
mean  translational  kinetic  energy  per  unit  mass  of  a  particle  of 


species  a  relative  to  mass  velocity  of  mixture; 

V 

—a 

velocity  of  a  particle  of  species  a 

V 

mass  velocity  of  mixture;  pv  =  Z  p  v 

—  Q  a—a 

Xa 

ambipolar  diffusion  velocity;  v^  =  Xp  " 

V 

V 

mass  velocity  of  plasma;  p^v^  =  Pj^v^  -f 

Pe^e 

V 

—a 

mass  velocity  of  species  a;  v  =  (1  /n  ) 

”*a  a 

/f  dV 
/  a — ( 

at 


at 


w 
—  a 


w 

-ap 


diffusion  velocity  of  species  a  relative  to  mixture;  -  v 

diffusion  velocity  of  species  a  relative  to  plasma; 

w  =  V  -  V 
—  ap  —a  — p 

energy  per  unit  mass  of  a  particle  of  species  a,  excluding 
translational  energy 

mass  density  of  mixture;  P  =  P 


mass  density  of  plasma;  p  =  P^  P^ 


mass  density  of  species  a;  p^  =  ^a^a 


shear  tensor  of  plasma  t  =1-  +  I 

^  =p  -ip  -ep 


shear  tensor  of  mixture  I-  =  Z  X 

t  n 


at 
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NOMENCLATURE  (Continued) 


T 

=a 


shear  tensor  of  species  a  relative  to  mass  velocity  of  species; 


T  =  P  I  -  P  c  c 
-a  a— a— a 


T 

=ap 


shear  tensor  of  species  a  relative  to  mass  velocity  of  plasma; 


T  =p  I-pC  C  =p  I-pCC  -  pw  w 
=ap  ^ap=  a— ap— ap  cip“  a— a— a  a— ap—ap 


=  at 


shear  tensor  of  species  a  relative  to  mass  velocity  of  mixture; 


I«.  =  p4.l“pC.C  =p  I-p  C^  -  p  w  w 
-at  ^at=  a— at— at  ^at=  a— a— a  a— a— a 


u> 


mass  rate  of  production  per  unit  volume  of  species  a 
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